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CHAPTER  I 


INTRODUCTION 

Composites  are  becoming  more  widely  used  as  primary  structural 
materials,  hence  the  concern  for  predicting  reliability  is.  also  increasing 
As  fiber-reinforced  composite  materials  make  the  transition  from  labora¬ 
tory  test  coupons  and  non-critical  components  to  major  production  items, 
methods  for  predicting  the  fracture  and  strength  characteristics  are 
necessarily  drawing  more  attention.  The  question  of  attempting  to  under¬ 
stand  the  nature  of  fracture  and  damage  growth  in  laminated  composites 
can  be  approached  in  several  different  ways.  One  may  analyze  the  laminate 
as  being  composed  of  a  homogeneous  orthotropic  material  (or  of  homogeneous 
orthotropic  layers)  and  use  numerical  methods  such  as  finite  elements 
to  obtain  results  for  the  response  of  particular  structures  or  even  to 
investigate  local  effects  near  geometric  or  load  discontinuities.  The 
use  of  an  equivalent  homogeneous  material  and  of  numerical  methods 
to  investigate  local  effects  is  not  always  satisfactory  [1].  It  is  often 
the  multiphase  nature  of  a  composite  that  gives  it  its  unique  properties, 
especially  near  the  damaged  regions.  To  replace  the  actual  material 
by  an  equivalent  homogeneous  material  may  well  remove  one  of  the  most 
important  characteristics  of  the  material. 

Another  approach  is  to  use  empirical  models  (descriptive  models) 
which  generally  do  not  account  for  the  physics  of  the  particular  materials 
but  contain  parameters  that  are  adjusted  to  gain  results  agreeing  with 
experiments.  Such  descriptive  models  have  little  use  in  understanding 
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damage  growth  and  local  effects,  but  can  be  valuable  in  developing 
design  guidelines  for  structural  components. 

By  contrast,  the  present  study  is  concerned  with  developing  rela¬ 
tively  simple  mathematical  models  that  contain  the  important  physical 
and  geometric  properties  of  the  composite  (predictive  models)  such  that 
it  is  possible  to  investigate  the  fundamental  behavior  of  a  laminate 
in  terms  of  the  various  material  properties  of  each  lamina.  One  of 
the  main  difficulties  with  such  modeling  is  selecting  significant 
properties  to  have  a  model  complete  enough  to  be  reasonably  accurate. 

At  the  same  time  it  should  be  simple  to  solve.  Zweben  gives  an  excel¬ 
lent  discussion  of  this  "materials  modeling  approach"  in  [2]. 

Some  of  the  first  work  in  modeling  a  unidirectional  composite  in 
this  manner  was  done  by  Hedgepeth  [3] ,  where  no  damage  other  than  an 
initial  transverse  notch  was  considered.  The  study  was  extended  by 
Hedgepeth  and  Van  Dyke  [4]  for  the  special  case  of  one  broken  fiber 
with  longitudinal  splitting.  The  extension  to  more  than  one  broken  fibe 
with  longitudinal  damage  could  not  be  developed  conveniently  by 
influence  functions  as  in  [3]  and  [4]  because  the  broken  fiber  adjacent 
to  the  damage  region  was  not  typical  of  any  of  the  remaining  fibers. 
Goree  and  Gross  [5]  used  Fourier  transforms  to  modify  the  solution  to 
account  for  an  arbitrary  number  of  broken  fibers  as  well  as  for  longi¬ 
tudinal  matrix  damage  to  include  both  yielding  and  splitting  initiating 
at  the  notch  tip  between  the  last  broken  fiber  and  the  first  unbroken 
fiber.  Goree,  Dharani  and  Jones  [6]  added  constraint  layers  to  the  main 
lamina  to  account  for  either  a  misalignment  of  fibers  in  a  multi-ply 
unidirectional  laminate,  or  for  the  presence  of  angle  plies  which  give 
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support  to  the  unidirectional  ply.  The  constraint  layer  was  taken  as 
being  fully  bonded  to  the  unidirectional  lamina  at  all  times  in  the 
above  model. 

The  analysis  presented  here  is  an  attempt  to  extend  the  above  dis¬ 
cussed  model  to  include  the  effect  of  debonding  between  the  notched 
unidirectional  lamina  and  the  surface  constraint  plies.  The  laminate 
is  modeled  as  a  two  dimensional  region  of  a  unidirectional  lamina 
with  symmetrically  located  surface  constraint  layers  whose  fibers  make 
an  angle  theta  with  the  unidirectional  ply,  (Figure  1). 

In  the  vicinity  of  a  notch  in  a  laminate,  the  broken  fibers  exert 
longitudinal  shear  stresses  in  the  matrix  which  are  transferred  to  the 
nearest  unbroken  fibers.  The  shear-lag  assumption  [2,  7,  8]  is  used 
for  this  shear  transfer  between  fibers  in  the  unidirectional  lamina. 

An  equivalent  expression  is  used  for  the  corresponding  effect  of  the 
constraint  layers.  Since  the  shear  transfer  does  not  depend  on  the 
transverse  displacements,  it  uncouples  the  longitudinal  and  transverse 
equilibrium  equations.  Hence,  the  longitudinal  displacements  can  be 
calculated  without  solving  the  transverse  displacement  equations. 

The  fibers  support  all  of  the  longitudinal  stress  in  the  unidirec¬ 
tional  lamina  because  the  longitudinal  extensional  modulus  of  the 
fibers  is  assumed  to  be  much  larger  than  that  of  the  matrix.  Debonding 
between  plies  in  the  vicinity  of  the  crack  is  considered  to  be  of 
finite  width  and  extends  to  infinity  in  the  longitudinal  direction. 

The  extreme  fibers  of  the  debond  zone  are  assumed  to  be  attached  to  each 
other  across  the  debonded  region  by  a  spring  of  stiffness  depending  on 
the  constraint  layer  properties  and  on  the  width  of  the  debonded  zone. 
That  is,  the  layer  debonds  but  is  still  connected  to  either  side  of 


fiber 


debond ins 


Figure  1. 


Unidirectional  Lamina  with  Broken  Fibers,  Surface 
Constraint  Layers  and  Debonding. 
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the  region  and  carries  load  due  to  the  longitudinal  displacement  of 
the  extreme  fibers  of  the  debonded  zone.  This  can  be  more  clearly  seen 
by  referring  to  Figure  1. 

As  an  initial  investigation,  the  basic  mechanism  of  crack  growth 
is  limited  to  a  model  containing  broken  fibers  only.  Subsequently,  a 
model  is  next  developed  to  account  for  additional  longitudinal  damage 
parallel  to  fibers  in  the  monolayer,  (Figure  2).  Splitting  and  yield¬ 
ing  of  the  matrix  is  assumed  to  initiate  at  the  notch  tip  and  to  pro¬ 
gress  longitudinally  between  the  last  broken  fiber  and  the  first 
unbroken  fiber. 

There  are  three  different  zones  in  the  model. 

(i)  unidirectional  ply  with  bonded  constraint  layers,  fiber 
numbers  (0)  to  (N-l)  and  (M+l)  to  (°°). 

(ii)  unidirectional'  ply  with  debonded  constraint  layers,  fiber 
numbers  (N+l)  to  (M-l),  and 

( iii)  intermediate  fibers  (N)  and  (M>  of  the  following  ply. 

An  equilibrium  equation  is  written  for  each  fiber  using  the  basic 
stress-displacement  relations  given  by  Hooke's  law  and  the  shear- lag 
assumption.  A  description  of  the  solution  to  these  equations  is  given 
in  the  following  sections.  The  stresses  and  displacements  are  determined 
as  a  function  of  number  of  broken  fibers,  constraint  layer  parameters 
and  debonded  zone  size.  The  results  are  compared  to  the  corresponding 
fully  bonded  cases. 


X 


Figure  2 


debonded  zone 


Unidirectional  Lamina  with  Broken  Fibers,  Longitudinal 
Matrix  Damage,  Surface  Constraint  Layers  and  Debonding 
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CHAPTER  II 
FORMULATION 

Two  Dimensional  Shear-Lag  Model  with 
Broken  Fibers,  Surface  Constraint 
Layers  and  Debonding 

A  unidirectional  array  of  parallel  fibers  with  surface  constraint 
layers,  debonding  and  an  arbitrary  number  of  broken  fibers  is  shown  in 
Figure  1.  Debonding  is  assumed  to  exist  from  the  last  broken  fiber  (L) 
to  an  unbroken  fiber  (M)  of  the  unidirectional  lamina.  The  constraint 
layers  are  intended  to  represent  adjacent  layers  of  a  unidirectional 
lamina.  They  are  assumed  to  be  placed  symmetrically  about  the  uni¬ 
directional  lamina  to  give  a  laminate  with  no  bending.  The  broken 
fibers  are  assumed  to  occur  along  the  x-axis  and,  since  the  loading 
is  symmetric,  only  the  first  quadrant  of  the  laminate  is  considered 
in  the  analysis.  The  basis  analysis  and  assumptions  are  the  same  as 
in  [6],  however,  in  order  to  clearly  indicate  the  modifications  needed 
to  account  for  surface  debonding,  it  is  necessary  to  repeat  the  basic 
formulation . 

The  fibers  are  taken  to  be  of  much  greater  stiffness  and  strength 

than  the  matrix  and  the  longitudinal  load  is  therefore  assumed  to  be 

carried  by  fibers  only.  Load  is  transferred  between  fibers  by  shear 

stresses  as  given  by  the  classical  shear-lag  assumption.  The  axial 

fiber  stress  a  (y)  and  the  matrix  shear  stresses,  t  (y)  and  TT(y), 
n  n  n 

are  then  given  by  the  simple  relations: 
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G T 

XA(y)  =  [vn(y)  "Vn-l(y)^  (3) 

where 

vn(y)  =  axial  displacement  of  the  fiber  (n)  at  the  location  (y) , 

Ep  =  Young's  modulus  of  the  fiber, 
t  =  thickness  of  the  unidirectional  ply, 
t'  =  thickness  of  the  constraint  plies. 

The  stiffnesses  G^/h  and  G^/h'  must  account  for  interaction 
between  fibers  [5,  7,  8].  and  G^  are  typically  not  the  shear  moduli 

for  the  "neat"  matrix  nor  are  h  and  h'  necessarily  fiber  center-line 
distances.  The  ratios  G^/h  and  G^/h1  are  equivalent  stiffnesses  and 
are  assumed  to  be  material  constants  and  depend  only  on  the  fiber  and 
matrix  properties,  the  fiber  volume  fraction,  orientation  of  plies, 
and  not  on  the  size  of  the  damage  region.  Only  for  large  spacing  can 
Gm  and  h  be  expected  to  approach  the  "neat"  matrix  and  center-line 
values . 

By  the  virtue  of  the  shear-lag  assumption  the  longitudinal  and 
transverse  equilibrium  equations  become  uncoupled  and  the  longitudinal 
displacement  and  stress  in  the  fibers  as  well  as  the  matrix  shear 
stress  can  be  obtained  without  solving  the  transverse  equilibrium 
equations.  Therefore,  only  the  equilibrium  equations  in  the  longi¬ 
tudinal  direction  will  be  considered  in  the  following  discussion. 
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As  pointed  out  earlier  and  in  [5],  G^/h  and  G^/h '  are  to  be  de¬ 
termine  experimentally.  For  example  in  [9]  and  [10],  it  is  shown  that 
the  shear  stress  becomes  larger  as  the  fiber  spacing  decreases,  that 
is ,  (0(1/ /d) )  for  rigid  fibers  where  Td!  is  the  minimum  distance  be¬ 
tween  the  fibers.  Local  failure  may  occur  at  critical  points  through 
the  thickness  in  advance  of  laminate  splitting  which  would  give  an 
apparent  shear  stiffness  considerably  different  from  that  of  the  matrix 
alone. 

The  debonded  fibers,  (N)  and  (M)  are  considered  to  be  connected 
by  springs  due  to  the  presence  of  the  angle-plies  of  the  constraint 
layers.  The  springs  are  assumed  to  have  a  linear  force-displacement 
relation  and  the  stiffness  (k)  per  unit  area  for  a  particular  laminate 
to  decrease  proportional  to  the  length  of  the  spring. 

With  reference  to  the  free-body  diagrams,  Figures  3  through  5,  of 
the  elements  for  different  ranges  of  fibers,  the  equilibrium  equations 
are 


K  d°V  ,  t' 

—  ,  i  —  T  )  +  (t,|l1-t'  )  V 

t  dy  n+1  n  |  n+1  n  t 


for  n  =  0,  1,  2,...,N-1,  M+l,..., 

K  dcV 

T  IF  +  <T|„+1-T|n)  ■  0 


for  n  =  N+1,  N+2,...,M-2,  M-l, 


where 

A  =  area  of  fiber. 

r 

For  fiber  N 


(4) 


(5) 


NOTE:  Diagram  shown  as  above  for  clarity  only.  In  the  actual  laminate 

the  constraint  layers  are  symmetric  with  a  layer  of  thickness 
t'/2  placed  on  each  side  of  the  unidirectional  ply. 


Figure  4.  Free-body  Diagram  for  a  Typical  Element  of  the  Fully  Bonded 
Zone  (Fibers  No.  (0)  to  (N-l)  and  (M+l)  to  (°°))- 


}  C  (y jOy)  1N 


T(y-J-Oy  ) 


,(y)lK 


^  / 

■*  / 

/  «EJ=> 


i  /^_jZl(y),w  |  | 

JZ_ 1 _ |  J4r 


/  ^(y)iK  v  — |- 


(y )  1  w 


NOTE s  Diagram  shown  as  above  for  clarity  only.  In  the  actual  laminate 
the  constraint  layers  are  symmetric  with  a  layer  of  thickness 
t ' / 2  placed  on  each  side  of  the  unidirectional  ply. 


Figure  5.  Free-body  Diagram  for  a  Typical  Element  of  the  Intermediate 
Zone  (Fiber  No.  (N)). 
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For  fiber  M 

t  dy 


+  (t 


M+l 


-  T 


M' 


)  - 


kt’ 


(v  —  V  )  +  T  * 
^  M  N; 


td 

M+l  t 


=  0. 


(7) 


Using  the  stress-displacement  relations  (1) ,  (2)  and  (3) ,  in  the 
above  equilibrium  equations,  the  following  set  of  difference-differen¬ 
tial  equations  is  obtained: 


A  Eh  d  v 

— — - —  +  (1  +  Cn)  (v  -  -  2v  +  v  )  =  0 

G  ,t  ,  2  v  R'  n+1  n  n-l/ 

M  dy 

for  n  =  0,  1,  2,...,N-1,  M+l,..., 

A_,E_h  d2v 

— - — - h  (V  ,  -  2v  +  V  ,  )  =  0 

G„,t  ,2  n+1  n  n-r 

M  dy 

for  n  =  N+1,  N+2 , . . . ,M-2 ,  M-l. 

For  fiber  N 


(8) 


(9) 


G  t  2  +  ^VN+l"  2vN  +  VN-1^  +  CR2('VM  VlP  CR^VN  VN-1^  °' 
M  dv 


M  dy 
For  fiber  M 


G  t  2  +  (VM+1  2VM  +  VM-1)  CR2('VM  VN)+CR('VM+1  VM^  °'  (11) 

M  dy 


where 

c  = _ i _ 

R2  (GM/h) 


C 


R 


(G’/h,)  t, 
(GM/h)  t  ‘ 


(12) 


The  constraint  layer  provides  additional  longitudinal  stiffness  to 
the  unidirectional  ply,  the  effect  being  given  by  the  constant  CR.  The 
debonding  effect  is  represented  by  the  second  constraint  ratio  * 
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To  match  the  differential  equations  for  the  bonded  case  [9], 
must  reduce  to  CR  when  M-N=l.  The  only  varying  parameter  for 
in  a  particular  laminate  is  the  width  of  the  debond  zone,  (M-N). 
Since  C  £  represents  a  linear  spring  and  loses  its  stiffness  (k)  per 
unit  area  proportional  to  its  length  then  C ^  * 


Hence, 


R2  M-N 


Noting  the  coefficient  of  the  second  derivative  in  Equations 
through  (12) ,  the  following  changes  of  variables  as  suggested  in 
are  made. 


(13) 

(8) 

[6] 


y 


n* 


dv 

n 

a  =  a  a  =  -z — * 

n  00  n  F  ay 


v 

n 


N 


Aph 


w  v 


where 


=  applied  remote  stress, 
and  q,  and  V^(n)  are  non-dimensional. 

Algebraic  manipulation  of  Equations  (3)  and  (14)  gives 


dV 

n 


(14) 


(15) 


t  =  a 
n  c 


gmaf 


'  AjEjTt  (vn  ~  Vn-T‘ 


(16) 


The  resulting  equations  in  non-dimensional  form  are 


d2v 


dq 


(1  +  CR>  (Vi-^n+Vd  = 0 


(17) 


for  n  =  0,  1,  2,...,N-1,  M+l,..., 


d2v 


n 


dn 


2  +  <Vl-2Vn+Vl>  '  ° 


for  n  =  N+l ,  N+2,...,M-2,  M-l. 


For  fiber  N 


d2v. 


dn 


T+  <Vi-2vb  +  Vi)  +  ce2<vm-V-cr(vk-vn-i)  ■  °* 


For  Fiber  M 


d2v 


M 


dn 


2  +  (Vl-2VM  +  Vl)-CE2<VM-V-CK(Vl-V  =  °' 


These  equations  can  be  written  as  follows,  where  the  left-hand 
side  is  the  same  in  each  equation. 


d2v 


f+  U  +  CR)(Vl-2Vn  +  Vn.1>  -  o 


dn 


for  n  =  0 ,  1,  M+l,..., 

T+  <1  +  CRHVn+l-2Vn+Vn-l)  '  VVfVVl’ 


d2V 


dn 


for  n  =  N+l,  N+2,...>M-2,  M-l. 


For  fiber  N 


d2V 


dn 


f  +  a  ■ +  V  C  Vl  -  2Vh  +  VN_J_)  -  Cr(Vn+1  -  V  -  CR2(VM  -  V 


] 


For  fiber  M 


d2v. 


M 


dn" 


+  (l+CR)(Vtffl-2VM+vM.1)  -CR(VM-VH_l)  +  CR2(VM- 


These  difference— differential  equations  may  be  reduced  to 
tial  equations  as  in  [6]  by  introducing  a  new  function  V(n,9) 
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(18) 


(19) 

(20) 


(21) 

(22) 


(,)•  <23> 

V-  (24) 

different 

defined 


as 
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v  (n) 

V (n , 0 )  =  — - - 1-  E  V  (n)  cos  (n0)  (25) 

n=l  n 

from  which 

2  7r- 

v  (h)  =  —  /  V(n,0)  COS  (n0)d6.  (26) 

LI  (I  _ 


Making  use  of  orthogonality  of  the  circular  functions,  Equations 
(21)  through  (24)  are  then  written  as  one  equation,  valid  for  all  n  and 
n  as 


2  ir  fd2V  -'l 

~  f  <  2  “  2(1  +  CR)  [1  -  cos  (6)  ]  V  l  cos  (n0)  d0 

0  ^dn  J 

2  17  f  M  1 

-  “  “  /  <  L  G  (n)  cos  (£0)  cos  (n0)  )  d6, 

0  (.£=N  *  J 

where 

-  ■VVr!,.«H) 

for  £  =  N+l, . . . ,M-1 , 

GN(n)  =  -  cr(vn+1  -  vN)  +  cr2(vm~  vn)  5 
GM^n)  =  CR^VM"  VM-1^  ”  CR2^VM“  VEp  s 

G0(n)  -  -  |  I  Vvi  “  V  ‘  CR2<V  V  >• 

The  equation  is  of  the  form 


(27) 


(28) 


(29) 

(30) 

(31) 


—  /  F(n,0)  cos  (n0)d0  =  0  for  all  n  and  n 
0 


(32) 


and  as  F(ri,0)  is  even  valued  in  0 ,  if  the  integral  is  to  vanish  for 
all  n,  the  function  F(r),0)  then  must  be  zero.  The  single  equation 
specifying  V(n,0)  is  then 


2-  M 

d  V  2- 

— -  _  S  V  =  -  Z  G.(n)  Cos  (£0) 
dr)  £=N 


(33) 


<5  =  2(1+ CR)  [1-  cos  (0)]. 


where 
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The  solution  to  the  problem  of  vanishing  stresses  and  displace¬ 
ments  at  infinity  and  uniform  axial  compression  on  the  crack  surface 
will  now  be  sought.  The  complete  solution  will  be  obtained  by  adding 
the  results  corresponding  to  uniform  axial  stress  and  no  broken 
fibers  to  the  solution  of  the  following  problem.  The  appropriate 
boundary  conditions  are: 

dv  (n) 

v  (r0  =  0  ,  - j -  =  0  as  n  + 00  for  all  fibers,  (34) 

n  dp 

dV  (n) 

— - -  =  a  (n)  =  -  1  at  n =  0  for  all  broken  fibers,  (35) 

dr)  n 


Y  (*)  =0  at  n  =  0  for  all  unbroken  fibers.  (36) 

nv 

The  complete  solution  to  Equation  (24)  ,  satisfying  vanishing  stresses 

and  displacement  at  infinity  is  given  by 

oo  M 

V(n,e)  =  A(0)e_(6n)  +  /  4  sinh(6(n-t))  2  G  (t)  cos  (£0) dt  ,  (37) 

n  £=n 

where  the  function  A(0)  is  yet  unknown.  The  remaining  two  boundary 
conditions  give 


dV  (0) 
n 

dn 


2  11 
-  f 

11  0  «- 


"f 

0  V 


M 


-5A(0)+ /  <cosh(St)  Z  G  (t)  cos  (£0)  dt }  cos(n9)|d0 

£=N 


t|  cos(n0)  ' 


(38) 


for  all  broken  fibers,  and 


V  (0)  =  -  /  Ta(0)  -T  f  /sinh(6t)  Z  G„  (t) cos  (£0)dt)  cos  (n9)l d0 

n  ^  o  L  °  0  \  £=N  J  -* 


=  0 

for  all  unbroken  fibers. 

Equation  (37)  can  be  solved  exactly  by  taking 

,  oo  M  L 

A(0)  -  y  /  sinh(6t)  Z  G5  (t) cos (£9)dt  =  Z 
6  0  £=N  m=0 


B 

m 


cos (m0)  . 


(39) 


(40) 


Eliminating  A(0)  from  Equations  (38)  and  (40),  the  stress  boundary 
condition  (35)  reduces  to 


2  Tf  L  03  M 

■T  ’  6  E  Bmcos(m  6)+/  e  V  '  Z  Gf (t) cos (£0) dt  ^  cos(n0)d0  (41) 


71* 


0 


£=N 


0  v  m=0 
=  -  1. 

From  Equations  (37)  and  (40)  ,  A(9)  can  be  eliminated  to  obtain 
V(r)50)  in  terms  of  B^  and  G0 (t)  .  Recalling  the  relation  between 


V(t]  ?  6)  and  V^(q)  an  expression  can  be  obtained  for  longitudinal  fiber 


displacements  as 


V  (q)  -  —  f  fe  Z  B  cos(m0) 

Ti  7T  i  m  s  ' 


L 


n 


m=0 


m 


1  °°  D  ( 5  t  ^ 

+  —  f  -  ~ ?  E  G. (t) cos (£0) dt  >  cos(nO)d0  (42) 

0  £=N  *  J 


where 


D(6,q,t)  =  e 


-6 1  n-t |  -6(n+t) 


-  e 


(43) 


The  longitudinal  fiber  stress  is  obtained  by  differentiating 

Equation  (42)  with  respect  to  q  and  is 

L 
Z 

0  m=0 

<5 1  n-t 


dV  (q)  tt 

=  -V  =  f 


-  —  /  (  -6  Z  B  cos(m0)e 
■k  \  m  x  7 


-(5q) 


+  / 


0  L 


-6(q+t) 


L 


E  G.  ( t)  cos  (£0) )  cos  (n0)d0, 
£=N  j 


where 


(44) 


p  =  1  for  t  <_  q  , 
p  = -  1  for  t  >  q  o 

Equations  (28),  (29),  (30),  (41)  and  (42)  can  be  solved  for  the  unknowns 


and  G. ( t) . 
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Limit  Case  of  an  Infinitely  Wide 
Debonded  Zone 

This  is  an  extension  of  the  model  developed  previously ,  such  uhat  the 
debonded  zone  is  now  assumed  to  extend  to  infinity.  Since  by  Equation  (12) 


tf 


JR2  (GM/h)  t 


and 


k  a 


r— ?  then  C  9  ■+  0 
M-N  R  2 


for  a  debond  zone  of  infinite  width. 

Physically,  the  spring  between  fibers  at  the  extremities  of  the  de¬ 
bonded  zone  has  no  stiffness  as  it  has  an  infinite  length.  Also,  for 
fibers  (M)  and  (M-l)  far  from  the  crack  tip  ~  Vm-1#  Then  usinB 

CR2  =  0  and  VM  =  VM-1 


Equations  (21),  (22),  (23),  (24)  reduce  to 
d2V 

~f+  <1+CR>  <Vl-2Vn  +  Vl>  -  0 

dn 

for  n  =  0,  1,  2,...,N-1,  M+l,...,  and 


(45) 


d2V 

n 


dn 


0  +  (V  ,  -  -  2V  +  V  -, )  =  0 
2  v  n+1  n  n-r 


(46) 


for  n  =  N+1,  N+2,...,M-2,  M-l. 
For  fiber  N 


d2v. 


N 


dp 


2^TvT  +  )  C-p  (VMa_1  V^)  , 


N  N-l 


Rv  N+1  N' 


(47) 


Following  the  same  technique  as  before,  the  single  differential 

equation  to  be  solved  is  then 

2“  N 

-  62  V  =  -  2  H  (n)cos(£6) 

dn  £=0 


(48) 
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where 


'  Wi-2vi+Vi) 


(49) 


For  fiber  (0) 


H0<">  *  CR(V0-V- 


(50) 


The  boundary  conditions  (34) 9  (35) ,  (36)  yield 


9  71  r  L  03  -(xt)  N  ^ 

—  /  <-S  £  B  cos(m9)  +  /  e  £  Hn  (t)  cos  (£0)  dt  )  cos(n0)d0  (51) 

"oIfo  “  0  *=0  1  J 


m 


=  -  1 

for  £  =  1,  25  3300e3N. 

Equations  (49) 9  (50) ?  (51) 9  (52)  can  be  solved  for  the  unknowns  B. 
and  H^(t) »  The  longitudinal  displacement  is  given  by 

V  (n)  =  -  /  fe~<'<5n')  Z  B  cos  (m0) 

n  77  n  1  n  m 

0  V  m“0 

1  ^  N 

+  y  £  H  (t)cos(£0)dt  )  cos(n0)d0,  (52) 

Z  0  °  £=0 


while  the  longitudinal  fiber  stress  is  given  by 


a  (ri)  =  —  /  f-S  £  Be  cos  (m0) 

n  7i r  n  1  n  m 

0  V  m=0 

TT  /  "6(n+t)  -6  I  n-t  I  \  N 

+;H — 


£  H  (t)cos(£0);  cos(n0)d0, 

'  £=0  *  J 

(53) 


where 


p  =  1  for  t  <_  T)  , 
p  =  -  1  for  t  >  n  . 
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Two  Dimensional  Shear-Lag  Model  with 
Broken  Fibers,  Longitudinal  Matrix 
Splitting  and  Yielding,  Surface 
Constraint  Layers  and  Debonding 

The  solution  developed  in  Section  I  will  now  be  extended  to  in¬ 
clude  longitudinal  splitting  and  yielding  of  the  matrix  as  shown  in 
Figure  2.  All  the  previous  assumptions  are  assumed  valid  and  it  is 
only  necessary  to  account  for  additional  damage.  It  is  assumed  that 
splitting  and  yielding  of  the  matrix  initiates  at  the  notch  tip  and 
progresses  longitudinally  between  the  last  broken  fiber  and  the  first 
unbroken  fiber.  The  matrix  is  assumed  to  be  elastic-perfectly  plastic. 
The  last  broken  fiber  is  considered  to  be  the  first  debonded  fiber. 

All  the  equations  remain  the  same  as  (4) ,  (5) ,  (6)  and  (7)  for  all 
fibers  except  for  fibers  (L)  and  (L+l).  The  equilibrium  equation  for 
fiber  (N)  is 


af  % 


N 


t  dy 
when  y  <_  , 

for  fiber  (N +1) 


-  To  <  y 


*2  > 


-  T 


.  kt T  ,  '  T 

N+  —  (vM-V  -T 


tj 
N  t 


=  0 


(54) 


S’ 


daT 


N+l 


dy 


+  T 


N+2  +  T0 


y  -  >  =  0 


(55) 


when  y  j<  . 

Tq  =  matrix  yield  stress, 

=  length  of  longitudinal  matrix  damage  at  the  crack  tip, 

=  length  of  longitudinal  matrix  split  at  the  crack  tip. 

The  above  equations  on  introduction  of  the  stress-displacement  rela¬ 
tions  (1),  (2),  (3),  become, 
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for  fiber  (N)  ,  when  y  <_  £ 


AFEF  d  VN  _  „  „  ^  GM  .  ^  .  kt'  , Tr  w  s 

t  2  T0<y  ~2  h  VN  VN-1  t  VM  V 

dy  r  t 

t ' 

-  T  <VN"VN-l)  T  "  °’ 


(56) 


and  for  fiber  (N  +  1)  ,  when  y  <_  £ 


1 


d  VN+1  GM  .  s  „  „  x  _  n 

— - ~ +  1T  (vn+2-vn+i)+to  <  y-^2  "  -  °- 

dy 


(57) 


Equations  (56)  and  (57)  are  re-written  as 


2 

AFEFh  d  VN 

IT  7Tt(1  +  V  <vn+i  "  2vn  +  vN-i> 

M  ay 


=  +  (1  +  CR)(VN+1“VN)-CR2(VM-VN)+T0  "  y'A2  > 


(58) 


AFEFh  d  VN+1  ,  .  w  „  s 

— Y~  +  (1  +  V  (vN+2  -  2vN+1  “  V 

M  dy 

=  -  (1  +  CR)  (vN+1  -  vN)  +  Cr(vn+2  -  vN+1)  -  x0  <  y  -  £2  >  .  (59) 


iW 


1  y 


a.  = 


2  =_N 


EFAFh 


G  t 
M 


(60) 


Vf 

T0  “  T0  ~\EFht 


The  resulting  overall  non-dimensional  equilibrium  equations  are 
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d2V. 


N 


dn 


+  <1  +  CR>  <Vl-2Vn  +  Vl>  '  0 


(61) 


for  n  =  0,  1,  2,...,N-lj  M+l , . . . . 
For  fiber  (N)  ,  when  r)  £  a 


d2V, 


dn 


#+  a+V(Vi-2vN+W 


-  <1  +  CR>  ‘VN+l-V-CR2<VM-V+'t0<r’- 


(62) 


For  fiber  (N) ,  when  n  >  a 


d2V. 


dn 


Y+  (1+Cr)(Vi-2?N  +  V»-1)-CR(V»+1-V-CR2<VM-V-  (63) 


For  fiber  (N+l)  ,  when  n  £.  a 


d2V 


dn 


2  +  O  +  CR) (VN+2  2vn+1  vn) 


=  -  a  +  V  (VN+1  -  VN)  +  Cr(Vn+2  -  VN+1)  -  r0  <  n  -  6  >.  (64) 


For  n  =  N+2,  N+3,...,M-1,  and  N+l,  when  n  >  a 


d2V 


n 


dn 


+  (l  +  CR)(Vn+1-2Vn  +  Vn-i)  -  CjtV^^  +  V^). 


(65) 


The  differential  equation  to  be  solved  is  the  same  as  Equation  (33) 
while  G  (n)  is  given  by 

V>  -  -tt+V(Vi-VtcJ!(v>-v()-;0  <n-  6>  for  me 


-  Wl- V  +  CR2<VM- V 


for  n  >  ot ,  (66) 


GN+i(n)  =  O  +  <  n  -  6  > 


R^VN+1  Ny  N+2  N+ly  '  l0 


S(Tl)  =  Wi-^  +  VF 

for  H  =  N+2,...,M-1  and  N+l  for  n  >  a. 


for  n  5_o i  ,  (67) 

(68) 
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Since  G  (n)  should  match  at  n  =  a  for  £  =  N+l  and  N, 


Va)  "  VN+l(a)  =  V 


(69) 


Equations  (41) 3  (42)  ,  (66) ,  (67) ?  (68)  and  (69)  can  be  solved  for  the 
unknowns  and  G^(t) „  The  expressions  for  displacements  and  stresses 
remain  as  given  by  Equations  (42)  and  (44) . 
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CHAPTER  III 
SOLUTION  TECHNIQUE 


In  all  the  problems  discussed  so  far,  the  solution  reduces  to  a 
series  equation  coupled  with  one  or  more  Fredholm  integral  equations 
of  the  second  kind.  Since  there  is  no  exact  form  of  solution  avail- 
able  to  solve  such  a  system  of  equations,  a  computer  program  is  de¬ 
veloped  by  modifying  a  numerical  procedure  given  in  [12] .  The  tech¬ 
nique  makes  use  of  a  method  by  Riez  [13]  to  solve  a  linear  Fredholm 
integral  equation  of  the  second  kind  defined  within  a  semi-infinite 
interval  of  integration.  A  given  integral  over  a  semi-infinite  in¬ 
terval  may  be  approximated  by  the  Gauss-Laguerre  quadrature  rule  as 
co  k  -x^ 

/  f(x)dx  =  E  w.e  f(x.)?  (70) 

0  i=l  1  1 

where  x.  is  the  i^1  zero  of  the  Laguerre  polynomial,  L,(x),  and  w. 

l  K  i 

is  the  corresponding  weight  function  given  by 


w. 

i 


x. 

l 


[(k+l)Lk+1(x.)]2 


(71) 


The  Laguerre  polynomial  L^(x)  is  given  by 


(72) 


Since  the  form  of  the  equations  for  each  solution  is  the  same  for 
the  three  particular  cases  discussed  above,  the  development  and  appli¬ 
cation  of  the  numerical  procedure  can  be  demonstrated  without  any  loss 
of  generality  by  taking  equations  corresponding  to  one  of  the  solutions 
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Consider  the  problem  of  an  infinite  debonded  zone  with  a  con¬ 
straint  layer  given  by  the  equations 


*¥">  -  Wi  -  2V» +  Vi> 


(49) 


cor  %  =  1,  2. 


,N-1,  N , 


H0(n>  -  CR(VrV0); 


(50) 


2  71  (  ^ 

—  f  I  -6  E  B  cos (m0) 

TT  n  |  n  m 

0  l  m=0 


0°  ,  ^  v  N  'S 

4-  /  e  0t  E  H«  (t)  cos  (£0)  dt  \  cos  (n0)d0 
0  £=0  ) 


-  1. 


(51) 


:or  n  =  03  1,...,L-1,  L9 


V  (q)  =  —  /  (e  E  B  cos(m0) 

n  77  n  1  n  m 

0  h  m=0 


N 


+  —  /  E  H  (t)  cos(£0)dt  }  cos(n0)d6„  (52) 

2  0  6  A=0  £  J 


Substituting  the  expression  for  V^(q)  in  (49)  and  (50)  yields 


Vn>  -  cR 


?  if  -('fin') 

-  /  e  E  B  cos  (m6) 

17  n  o  m 

0  m=0 


X  {cos  (  (£  +  1)  6)  -  2  cos(Jie)  +  cos  ( ($,  -  1)  0)}d0 

+  -//  D(~lS’T1,t')  E  H  (t)cos(nO) 

T  0  0  6  n=0  n 

X  {cos((£.+ 1)6)  -  2  cos(£6)  +  cos((£  -  1)0)  }d6dt 
L 


,(73) 


H0(n)  =  CR 


“  /  e  E  Bmcos(m0)  {  cos  (0)  -  l}d0 

'  0  m=0 


N 


+  —  /  /  ^  E  H  (t)cos(n0)  {cos(0)  -  l}dtd0 

71  0  0  n=0  n 


(74) 
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The  integral  over  the  spatial  variable  't'  in  the  Equations  (51),  (73) 
and  (74)  can  now  be  replaced  by  the  series  in  Equation  (70)  to  yield 


L  2C  it 
E  — ^  / 
m=0  *  0  L 


•  -(<$n.) 

e  cos  (m0)  (  cos(  (£  +  1)  9)  -  2  cos  (5-0) 


4-  cos  ((&-!)  0)  d0 


+ 


6..-—  E  (k*  1(n.,t  )  -  2Kj(n±,t .) 
ij  it  j=1  ^  1  1  3  £  1  3 


t . 


+  \  (VVj  6  J  Wt.  VV 


C-n  N  k  r  ^ 

—  EE  {  K*+i(n,,t,) -2K  (n,  ,t.) 
•’  n=0  j-ll  n  1  J  "  1  3 


for  i  =  1,2,...,N-1,  N, 


+  K^1(p.,t.)j  e  3  wt_  Hn(t.) 


=  0,  (75) 


l  2C„  IT  r  -(5n.) 


E  —  / 
m=0  77  0 


1  cos  (m0)  {  cos  (0)  -  1}B 


m 


+ 


+ 


d0 


t . 


}  WVV'^VV}  e  J  Wt.  W 


13  7T 


3=1 


3 


C_  N  k  [  ,  0  It. 

J  E  E  (Kn(„.;t  )  -  Ka(n±,t  )|  e  »  V^) 

n=l  j=l  k  J  3 


0,  (76) 


and 


L  „  it 

Z  —  /  [-  6 cos(m0)cos(n0)B  ]d0 

TT  U1 

m=0  0 

N  k  0  irr-(5t,)  >|  t. 

+  E  E  —  /  /e  cos  (5-0)  cos  (n0)  }  e  w 

*=0  j.l’ol  J  3 

for  n  =  0,1,...»L-1,  L, 
where 

7T 

Km(ri  t.)  =  f  D(5,n.  »t  .)cos(m0)cos(n0)d0, 

n  1  3  0  1  J 


G£(t.)d0=-1 


(77) 


(78) 


28 


The  solution  therefore  reduces  to  solving  a  system  of  linear 

algebraic  equations  for  the  constants  B  and  for  explicit  values  of 

m  r 

the  functions  H^(t)  at  specified  quadrature  points.  The  above  system 
of  equations  is  solved  by  the  method  of  Gauss-elimination  with  partial 
pivoting. 

The  above  procedure  is  followed  in  all  the  three  cases;  however , 
for  the  instance  of  longitudinal  matrix  damage,  an  additional  step  is 
required.  The  presence  of  longitudinal  damage  has  a  very  significant 
effect  on  displacement  in  the  damage  region  and  hence  on  G^(t) .  If 
the  semi-infinite  integrals  are  expressed  as  a  series  expansion  given 
by  (70),  the  number  of  quadrature  points  lying  in  the  damage  region  is 
not  sufficient  to  represent  the  displacement  function  accurately  in 
that  region.  For  example ,  out  of  35  quadrature  points,  only  five  lie 
between  zero  and  two ,  a  typical  value  for  alpha  (a).  So  in  order  to 
approximate  the  integral  more  accurately,  more  points  are  required  in 
the  region  from  zero  to  alpha  (a).  The  following  procedure  is  used. 

The  integral  is  split  as  sum  of  two  integrals 

00  Ot  oo  0t  03 

I  f(x)dx  =  f  f(x)dx  +  /  f(x)dx  =  /  f(x)dx  +  I  f(x  +  ct)dx.  (79) 

0  0  a  o  0 

The  finite  interval  integrals  can  be  approximated  by  Gaussian  integra¬ 
tion  [14]  as 

a  .  k* 

/  f(x)dx  =  -|  I  w*  f(y*),  (80) 

0  j=l  2  3 

where 


y*  =  q  +  l)x*. 


(81) 


x*  is  the  jth  zero  of  the  Legendre  polynomial  P  (x) , 
and  w*  is  the  corresponding  weight  function  given  by 
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[p;(xp]2. 


(82) 


The  Legendre  polynomial  is: 


,n 


P  (x)  =  — 
n  „n 


,2  1N.n 

(x  -  1)  . 


2  n!  dx 


(83) 


Physically  it  would  be  more  direct  to  specify  the  applied  stress, 
number  of  broken  fibers,  L,  and  determine  the  damage  zone  a  and  g  de¬ 
pending  on  given  yielding  and  splitting  conditions .  As  a  and  g  appear 
in  the  limits  of  the  integrals,  this  is  not  convenient  computationally 
and  it  is  easier  to  specify  the  number  of  broken  fibers,  N,  and  the 


damage  zone  parameters  a  and  g  and  to  compute  the  required  applied 
stress  a 

oo 


Ed! 
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CHAPTER  IV 
RESULTS 

Debonding  with  no  Longitudinal  Matrix 
Damage  in  the  Unidirectional  Ply 

ihe  effect  of  the  width  of  the  debonded  zone  was  of  particular 
significance  in  this  study.  Results  are  given  in  Figure  6  for  various 
numbers  of  broken  fibers  with  a  constant  constraint  ratio.  Debonding 
was  assumed  to  start  at  the  last  broken  fiber  and  extend  longitudinally 
to  infinity,  (Figure  2).  The  critical  fiber  is  defined  as  the  fiber 
xtfhich  has  the  maximum  stress.  The  stress  in  the  critical  crack-tip  fiber 
decreased  initially  for  a  small  debonded  zone,  but  subsequently  increased 
wii_h  an  increase  in  the  width  of  the  debonded  zone.  In  fact,  the  stress 
concentration  in  the  limit  case  of  an  infinitely  wide  debonded  zone, 
was  more  than  that  of  the  bonded  case,  (Table  I). 

One  result  of  particular  significance  is  that  the  maximum  decrease 
in  the  stress  in  the  critical  fiber  occurs  for  a  small  debonded  zone 
and  is  essentially  independent  of  the  initial  crack  length.  Figure  6 
shows  that  a  debonded  zone  width  of  two  fibers  spacings  results  in  the 
largest  decrease  in  the  maximum  fiber  stress  for  five,  seven,  nine  and 
twenty-one  broken  fibers. 

Debonding  acts  like  a  constraint  (CR2)  between  the  last  broken  fiber 
(N)  and  last  debonded  fiber  (M),  resulting  in  a  redistribution  of  stresses 
in  the  vicinity  of  the  crack  tip,  hence  decreasing  the  stress  in  the 
critical  fiber.  A  higher  constraint  ratio  resulted  in  a  larger  drop 
m  critical  stresses  as  shown  in  Figure  7  but  gave  higher  critical 
stresses  in  the  limit  case  Figure  8. 
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For  case  (ii),  unidirectional  lamina  with  no  constraint,  it 
was  found  in  [5]  that  once  the  split  forms  the  critically  stressed 
fiber  unloads  and  the  split  length  becomes  unbounded  under  a  five-to 
ten-percent  further  increase  in  applied  stress.  The  fracture  reduces 
to  an  unnotched  laminate  with  the  net-section  fracture  stress  being 
independent  of  the  initial  crack  length. 

For  cases  (iii)  and  (iv),  the  critically  stressed  fiber  does  not 
reduce  to  a  net  section  state  but  continues  to  carry  load  after  splitting 
with  increasing  applied  load.  However,  in  the  presence  of  debonding, 
the  maximum  fiber  stress  is  relieved. 

Hence,  in  terms  of  load  carrying  capacity,  the  worst  case  is  (i) 
where  the  monolayer  has  no  damage  other  than  an  initial  transverse 
notch,  while  the  best  case  is  (ii)  where  the  notched  monolayer  is  not 
constrained  but  has  longitudinal  matrix  damage.  Cases  (iii)  and  (iv) 
lie  between  the  above  models  where  debonding  allows  the  larger  load 
carrying  capacity. 

For  all  fully  bonded  constraint  layer  cases,  the  maximum 
fiber  stress  occurs  in  the  first  unbroken  fiber  at  the  end  of  the  split 
(y  =  )  for  no  (or  low)  constraint  ratio  [5]  and  at  the  notch 

tip  (y  =  0)  for  high  constraint  ratios  [10].  The  same  behavior  occurs 
for  debonded  cases,  but  in  case  of  high  constraint  ratios  and  high  values 
of  alpha  (a)  and  beta  ( 6 ) ,  the  maximum  fiber  stress  occurs  in  the  last 
debonded  fiber  at  y  =  0.  This  shows  that  under  the  above  conditions 
the  first  unbroken  fiber  is  highly  relieved  of  stresses  and  can  result 
in  discontinuous  damage  of  the  fibers.  The  following  table  quantita¬ 
tively  illustrates  the  above  behavior. 
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CHAPTER  V 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  main  aim  of  this  research  project  was  to  investigate  the  effects 
of  interlamina  debonding  between  notched  unidirectional  lamina  and  con¬ 
straint  plies.  The  significant  results  are  as  follows: 

(i)  The  maximum  stress  in  the  crack  tip  fiber  decreases  initially 
for  a  small  debonded  zone,  but  increases  as  the  width  of  the 
debonded  zone  is  increased. 

(ii)  The  largest  decrease  in  the  maximum  fiber  stress  occurs  for 

a  debonded  zone  width  of  two  to  three  fiber  spacings.  Further 
this  debond  zone  width  is  essentially  independent  of  crack 
length . 

(iii)  For  a  small  debond  zone,  the  maximum  fiber  stress  in  the  mono- 
layer  decreases  as  the  constraint  ratio  increases. 

(iv)  In  the  presence  of  longitudinal  matrix  crack- tip  damage  and 
constraint  layers,  the  maximum  fiber  stress  in  the  unidirec¬ 
tional  ply  is  further  relieved  if  debonding  takes  place. 

(v)  For  a  large  constraint  ratio  with  longitudinal  matrix  damage 
at  the  crack  tip,  debonding  reduces  the  stress  in  the  first 
unbroken  fiber  such  that  the  maximum  fiber  stress  can  occur 
at  the  end  of  the  debond  zone. 

The  following  modifications  are  recommended: 

Extend  damage  to  the  whole  laminate,  that  is,  account  for  damage 
in  the  constraint  plies. 

Represent  finite  longitudinal  debonded  zones  around  the  crack. 

Understand  the  concept  of  constraint  ratio  in  terms  of  geo¬ 
metric  and  physical  layer  properties.  Since  the  constraint 
ratio  is  not  defined  in  terms  of  layer  properties,  it  is  not 
possible  to  make  quantitative  predictions  about  laminate 
strength.  Approximate  best  fit  curves  [6]  and  experimental 
studies  [15]  can  help  in  determining  constraint  layer 
parameters  to  find  constraint  layer  ratios. 


(i) 

(ii) 

(iii) 
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16  Abstract 


This  study  considers  a  method  of  analysis  for  predicting  the  fracture 
behavior  of  a  notched,  unidirectional  lamina  in  the  presence  of  surface  constraint 
layers  with  debonding  between  the  unidirectional  ply  and  the  constraint  layers. 

Two  particular  cases  are  presented,  the  first  being  a  debonded  zone  of. finite  width 
with  no  longitudinal  damage  in  the  unidirectional  ply.  This  solution  is  then 
extended  to  include  longitudinal  matrix  yielding  and  splitting  in  the  unidi rectiona 

ply  at  the  crack  tip.  _  . 

The  analysis  is  based  on  a  materials  modeling  approach  using  the  classical 
shear-lag  assumption  to  describe  the  shear  transfer  between  fibers.  The  fracture 
behavior  of  the  laminate  is  studied  as  a  function  of  initial  crack  length,  con¬ 
straint  ratio  depending  on  the  relative  physical  and  geometric  properties  of  the 
constraint  plies  and  the  unidirectional  lamina  (not  quantitatively  predictable  at 
present),  and  width  of  the  debonded  zone.  The  results  indicate  that  debonding  can 
reduce  the  maximum  fiber  stress  at  the  crack  tip  on  the  order  of  ten  percent.  This 
effect  is  maximum  for  a  debond  width  of  two  or  three  fiber  spacings  and  is  indepen¬ 
dent  of  the  initial  crack  length.  As  the  debond  width  grows  beyond  this  point,  the 
maximum  stress  increases.  For  widths  of  about  ten  fiber  spacings  or  more,  the 
maximum  fiber  stress  is  larger  than  for  the  fully  bonded  case. 

In  the  presence  of  longitudinal  matrix  damage  the  same  general  behavior  is 
found;  however,  the  location  of  the  maximum  fiber  stress  is  quite  complex.  In  some 
cases  with  large  matrix  damage  and  a  high  constraint  ratio,  the  maximum  fiber 
stress  can  occur  at  the  end  of  the  debonded  zone  away  from  the  crack  tip. 
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